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1.  INTRODUCTION 


The  following  sections  present  work  completed  under  the  UNM/AFRL  research  grant  FA9453- 
11-1-0276  since  May  2012.  During  this  time  period,  the  majority  of  the  research  completed  was 
the  development  and  implementation  of  analog  physics  models  and  advanced  single-event 
physics  models  in  Geant4  [1],  A  Geant4  model  of  the  CEASE  detector  was  completed  that  will 
be  used  for  future  testing  and  simulations. 

2.  BACKGROUND 

In  Monte  Carlo  simulation  of  the  transport  of  charged  particles,  the  analog  and  advanced  single¬ 
event  physics  models  require  total  and  differential  interaction  cross-sections,  which  are  used  to 
sample  the  particle  distance  to  collision,  scattering  angle,  and  energy-loss.  One  example  of  the 
transport  equation  for  electrons  that  can  undergo  angular  deflection  and  energy-loss  interactions 
is  given  by 


Q  Vv/(£,D)  +  2^(£,D)=  f  2„(A,/i0)V/(A',Q')dD'  +  2^  (E+Q,Q)iy(E+Q,a)dQ.  (1) 

J4^  JQmm 

Solution  to  this  equation  using  the  Monte  Carlo  method  will  give  an  exact  solution  within 
statistics.  In  a  given  material,  distance  to  collision  is  sampled  from 

P(s)ds  =  X,  exp(-X,s)tfo.  (2) 

Angular  deflection  and  energy-loss  are  sampled  from  their  respective  cross-sections  that  are  both 
characterized  by  the  Rutherford  cross-section  in  their  most  simple  form.  The  Rutherford  cross- 
section  for  angular  deflection  is 


CjE) 

[l-ft+2 rKE)f  ' 


The  Rutherford  cross-section  for  energy-loss  is 


(3) 


^Q)=w  (4) 

where  C(E)  and  K(E)  are  material  constants.  There  are  several  cross-sections  that  are  more 
accurate  than  the  Rutherford  cross-sections  and  typically  used  in  practice  (like  the  Mott  cross- 
section  and  the  Moller  cross-section).  Regardless,  the  Rutherford  cross-sections  in  eq.  3  and  eq. 

4  illustrate  the  peakedness  of  the  cross-sections  near  forward  angles  and  small  energy  losses.  In 
addition  the  total  cross-sections  are  large.  Cross-sections  of  this  nature  make  Monte  Carlo 
simulations  extremely  computationally  inefficient  and  impractical,  as  detailed  in  the  statement  of 
work  [6],  The  advanced  single-event  physics  models,  referred  to  in  the  following  sections  as  the 
Generalized  Boltzmann  Fokker-Plank  method  (GBFP),  are  based  on  powerful  results  from  Lewis 
theory  [2]  and  higher-order  Fokker-Plank  expansions  [8]  [9],  To  summarize  the  connection  of 
the  theory  to  the  GBFP  method,  nearby  or  approximate  forms  of  the  analog  TE  are  formed  using 
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approximate  forms  of  the  cross-section.  The  approximate  cross-sections  should  preserve  a  finite 
number  of  angular  and  energy-loss  moments  of  the  analog  cross-section.  Increasing  the  number 
of  moments  preserved  results  in  an  increasingly  more  accurate  approximation.  This  approximate 
transport  equation  is  solved  using  existing  single-event  Monte  Carlo  simulation  techniques.  The 
moments  are  defined  by 


Afl  =  2xj'_Pl(MiZ,(r,E,/W/i  (5) 

and 

>',=Ie“eX('%E.0)<e  (6) 

where  Pj(fi)  is  the  Legendre  polynomial  of  degree  1  and  Qk is  a  power  moment.  The  cross- 
sections  used  in  the  GBFP  method  are  obtained  by  forming  a  system  of  equations  to  determine 
discrete  points  and  weights  or  function  parameters  such  that  a  finite  number  of  moments  are 
preserved.  For  example,  to  obtain  a  discrete  GBFP  cross-section  the  following  system  is  formed 

where  is  a  discrete  cross-section  given  by 


a  .(E,M)  =  'ZC‘^lMo-Ll  (8) 

n=i  2/r 

Inserting  this  expression  into  eq.  (7)  gives 

M,  =!>;<#„)«„<£).  (9) 

71=1 

The  system  in  eq.  (9)  is  a  2N  by  2N  system  where  theM,  ’s  are  assumed  to  be  known  and  the 
discrete  points  and  weights,  £„and  an,  are  determined  using  Gaussian  Quadrature  methods. 
Details  on  the  algorithm  to  obtain  a  discrete  cross-section  are  given  in  the  following  section.  The 
discrete  cross-sections  are  then  used  in  Monte  Carlo  simulation  to  obtain  dose  deposition  and 
angular  distribution  results  and  compared  to  the  exact  analog  solutions. 
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3.  METHODS,  PROCEDURES,  AND  ASSUMPTIONS 


3.1.  GBFP  Discrete  Cross-section  Generation 

Over  the  past  reporting  period,  a  method  for  obtaining  a  solution  to  eq.  (9)  was  implemented.  A 
solution  to  this  system  is  necessary  for  generating  angular  scattering  and  energy-loss  cross- 
section  data.  The  cross-section  data  are  then  used  to  sample  deflection  cosine  and  energy-loss. 
The  solution  to  the  system  in  eq.  (9)  is  similar  to  a  Gauss  quadrature  system  for  an  arbitrary 
weight  function.  Therefore,  the  method  used  to  obtain  a  discrete  cross-section  for  this  work  is  an 
application  of  several  Gauss  quadrature  algorithms.  One  procedure  for  obtaining  a  solution  to 
system  in  eq.  (9)  is  outlined  below 

1.  Use  the  Modified  Chebyshev  Algorithm  (MCA)  [3]  to  map  coefficients  of  classical 
orthogonal  polynomials  to  coefficients  of  polynomials  for  an  arbitrary  weight  function. 

2.  If  using  Radau  quadrature  (fixes  one  endpoint  at  1),  update  Jacobi  matrix  according  to 
Golub’s  description  [5], 

3.  Apply  Golub  and  Welsch  algorithm  [4]  to  obtain  points  and  weights  that  are  given  by  the 
eigenvalues  and  eigenvectors  of  the  Jacobi  matrix. 

The  MCA  is  used  to  map  the  recurrence  coefficients  ( a,/3)  of  a  known  classical  orthogonal 
polynomial,  ?r(x),  given  by 


XJ+i(X)  =  (X~  aj -  PjXj- 1(*)  ( 1  °) 

to  the  recurrence  coefficients  (a,b)  of  an  arbitrary  set  of  polynomials,  pt(x),  that  are  given  by 

Pm(x)  =  (x~ ai)Pi(x) ~ hiPi  iO) •  ( 1 1 ) 

The  polynomials,  p^x),  are  orthogonal  with  respect  to  some  arbitrary  weight  function,  w(x), 
and  the  modified  moments  of  this  weight  function  are 

Zj  =  £  )w(x)  (12) 

The  modified  moments  and  the  coefficients  of  the  classical  orthogonal  polynomials  are  used  as 
input  to  the  MCA.  The  algorithm  uses  mixed  moments, 

r  b 

°kj  =< Pk\Xi>=)m  Pk{x)nl{x)w{x)dx ,  (13) 

to  relate  the  coefficients  of  the  classical  orthogonal  polynomials  to  the  desired  orthogonal 
polynomials.  Once  the  coefficients  ( a,b )  are  known  the  Jacobi  matrix, 
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aJ 

0 

A 

aJ+ 1 

-Jbj+i 

0 

1 

<*j+ 2 

0 

.  0 

o 

\j^N  2  °W-1  \[^N- 1 

0  a/v7 


\ 


y 


(14) 


is  formed.  The  eigenvalues  of  the  Jacobi  matrix  correspond  to  the  discrete  points  and  the  first 
entry  of  the  eigenvectors  squared  corresponds  to  the  discrete  weights.  If  a  Radau  set  is  to  be 
used  the  matrix  given  in  eq.  (14)  must  be  modified  according  to  Golub’s  description  [4],  The 
following  table  gives  an  example  of  a  discrete  cross  section  for  1  MeV  electrons  in  Carbon. 

Table  1:  Four  Point  Discrete  Angular  Deflection  Cross  Section  Data 


Discrete  Deflection  Cosine  Point  ( ) 

Discrete  Deflection  Cosine  Weight  ( an ) 

-8.24443 1439774413e-01 

2.010871029870031 

1. 88547858 1349338e-01 

8.343040064585361 

5.609457938047900e-01 

5.17935159877821 le+01 

9. 82989761 00 14093e-01 

7.507247294868574e+03 

The  cross  section  given  in  table  1  should  preserve  up  to  2N-1  moments  of  the  analog  cross 
section  since  Radau  quadrature  can  integrate  polynomials  up  to  order  2N-1.  The  following  table 
gives  moments  of  the  analog  cross  section,  moments  calculated  using  the  discrete  cross  section, 
and  the  relative  difference  between  the  two  for  the  cross  section  in  table  1 . 

Table  2:  Moments  of  Four  Point  Discrete  Angular-Deflection  Cross-Section. 


Order  of  Moment 

Analog  Moment 

Discrete  Moment 

Relative  Difference 

0 

9.293256144065792 

9.293256144065797 

501.0744e-018 

1 

9.291615893044463 

9.291615893044472 

877.0350e-018 

2 

9.288794923713871 

9.288794923713887 

1.7546e-015 

3 

9.284946297381490 

9.284946297381518 

3.0091e-015 

4 

9.280171953381741 

9.280171953381785 

4.7669e-015 

5 

9.274548248431435 

9.274548248431500 

7.1547e-015 

6 

9.268136174132964 

9.268136174133055 

9.9230e-015 

7 

9.260986466855453 

9.261386787328156 

43.2265e-006 

Notice  that  the  relative  difference  in  the  moments  is  exact  to  machine  precision  until  the  order  7 
moment  is  reached.  At  this  point,  the  higher-order  moments  are  approximated  in  terms  of  the 
lower  order  moments. 
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The  purpose  of  this  method  is  to  extend  the  mean  free  path  of  the  particle  while  preserving 
moments  of  the  cross-section.  The  cross-section  given  in  table  1  began  as  an  order  five  cross- 
section  with  one  discrete  point  at  H  —  1 .  This  point  corresponds  to  a  non-interaction. 
Specifically,  it  is  the  result  of  a  collision  with  a  nucleus  that  does  not  change  the  direction  of  the 
particle.  Moreover,  the  probability  of  this  interaction  is  several  orders  of  magnitude  more 
probable  than  for  jj.  <  0.95.  The  purpose  of  obtaining  a  cross-section  with  a  point  at  fi  =  1  is 
because  this  point  can  be  removed  as  it  cancels  with  a  term  on  the  left-hand-side  of  the  TE.  Once 
the  point  and  weight  is  removed  the  remaining  weights  are  renormalized.  The  renormalized  total 
cross-section  is  substantially  smaller  than  the  analog  total  cross  section.  Hence,  the  mean  free 
path  of  a  GBFP  particle  is  extended.  The  process  of  removing  the  peaked  portion  of  the  cross 
section  in  a  consistent  manner  is  a  regularization  process. 

The  ratio  of  the  mean  free  path  of  the  discrete  cross-section  to  the  mean  free  path  of  the  analog 
cross-section  is  used  as  a  rough  estimate  of  the  speed-up  gained  by  the  GBFP  method.  The 
following  figure  gives  this  ratio  for  electrons  of  various  energies  in  carbon.  The  ratio  of  the 
mean  free  path  associated  with  angular  deflection  is  given. 


H 


Number  of  Discrete  Points  and  Weights 


Figure  1.  Ratio  of  Mean  Free  Path  of  GBFP  Cross-Section  to  the  Analog  Cross-Section. 

As  seen  in  fig.  1 ,  the  largest  gains  in  efficiency  are  at  higher  energies  where  the  cross-section  is 
substantially  more  peaked.  As  the  number  of  points  and  weights  increases,  the  ratio  decreases 
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because  as  more  points  are  added,  the  weight  of  the  point  near//  =  1  decreases  to  correctly 
preserve  moments.  Figure  1  also  indicates  that  there  is  a  clear  trade-off  between  accuracy  and 
speed-up.  However,  previous  studies  [7]  demonstrated  that  discrete  cross-sections  with  just  one 
angular  deflection  and  energy-loss  point  and  weight  are  accurate  for  integral  quantities  like  dose. 

3.2.  IMPLEMENTATION  OF  THE  GBFP  DISCRETE  MODEL  IN  GEANT4 

Geant4  is  a  toolkit  for  the  simulation  of  the  passage  of  particles  through  matter.  Its  areas  of 
application  include  high  energy  nuclear  and  accelerator  physics,  medical  physics,  and  space 
physics  [1].  A  general  use  of  Geant4  entails  developing  an  application  with  at  least  three 
mandatory  user-defined  classes  that  are  inherited  from  the  base  classes 
G4VDetectorConstruction,  G4VPhysicsList,  and,  G4VPrimaryGenerator Action.  Creation  of  a 
G4VDetectorConstruction  inheriting  classes  is  used  to  set  the  problem  geometry.  In  this 
inheriting  class,  the  dimensions  and  material  properties  for  the  model  are  set.  A  G4VPhysicsList 
inheriting  class  is  used  to  set  the  various  physics  processes  that  are  applicable  for  the  particles 
required  by  the  simulation.  In  this  inheriting  class,  the  processes  and  the  associated  models  are 
created  and  applied  to  the  appropriate  particles.  Finally,  creation  of  a 

G4VPrimaryGeneratorAction  inheriting  class  is  used  to  define  the  particle  source.  The  position, 
energy,  particle  type,  and  direction  are  set  in  this  class  (distributions  for  each  of  these  parameters 
can  be  used  as  well).  It  is  at  the  level  of  processes  and  models  where  the  GBFP  method  is 
implemented.  A  process  and  model  must  be  written  for  each  model  (Discrete,  Continuous, 
Hybrid)  and  for  each  interaction  (angular  deflection  and  energy-loss). 

The  implementation  of  the  GBFP  method  occurs  at  the  level  of  physics  processes  and  models. 
For  a  typical  user,  these  classes  are  not  modified.  The  user  selects  from  a  list  of  processes  and 
models  and  sets  which  physics  will  be  used  in  the  physics  list.  However,  creation  of  processes 
and  models  for  the  GBFP  method  will  require  additional  development.  Analog  and  discrete 
GBFP  processes  and  models  are  completed  to  date.  The  associated  classes  inherit  from 
G4VEmProcess  and  G4VEmModel.  The  process  class  controls  the  ordering  of  interactions.  For 
example,  a  process  can  be  applied  at  the  beginning  of  a  step,  along  a  step,  and  post-step.  The 
process  classes  require  an  associated  model.  The  models  handle  calculation  of  the  total  cross 
section  for  sampling  distance  to  collision  and  methods  for  sampling  quantities  like  angular 
deflection  and  energy-loss.  The  benefit  of  using  Geant4  is  gained  through  the  modularity  in  the 
code  package.  As  discussed,  the  implementation  of  the  GBFP  method  only  requires  a 
description  what  interactions  can  occur  and  how  to  simulate  these  interactions.  Beyond  that,  the 
code  handles  all  of  the  other  details  and  contains  additional  functionality  like  classes  for 
extracting  results  and  a  sophisticated  geometry  modeling  capability.  The  following  figure  shows 
a  three-dimensional  CEASE  model  captured  with  the  Geant4  visualization  tools.  With  the  ability 
to  easily  adapt  GBFP  models  and  model  sophisticated  geometries  for  testing  of  the  GBFP 
method,  Geant4  is  a  clear  choice. 
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Figure  2.  3-D  Rendering  of  Geant4  CEASE  Detector  Model. 
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4.  RESULTS  AND  DISCUSSION 


The  following  section  presents  results  from  initial  testing  of  the  GBFP  discrete  Model.  The  first 
several  figures  are  angular  distributions  obtained  using  the  GBFP  discrete  model  and  the  analog 
model.  Angular  distributions  were  calculated  for  1  MeV  and  100  keV  electrons  in  gold  slabs.  As 
seen  in  fig.  3,  the  discrete  model  cannot  resolve  angular  distributions  with  two  discrete  angles  in 
thin  slabs. 
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Figure  3.  Comparison  of  Angular  Distribution  for  1  MeV  Electrons  in  Gold. 

However,  increasing  the  number  of  angles  increases  the  accuracy  of  the  model  and  with  enough 
points  the  discrete  model  can  resolve  angular  distributions.  Figure  4  presents  a  comparison  of 
angular  distributions  for  1  MeV  electrons  in  gold  generated  using  a  32-point  discrete  model  and 
the  analog  model.  The  discrete  solution  is  in  good  agreement  with  the  analog  solution. 


Figure  4.  Comparison  of  Angular  Distribution  for  1  MeV  Electrons  in  Gold. 
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Adding  more  points  will  improve  the  accuracy  of  the  model.  This  indicates  that  the  algorithm 
for  generating  the  cross-section  data  and  the  implementation  of  the  model  was  done  correctly. 
Although  these  results  are  satisfactory,  there  are  a  couple  limitations  to  the  discrete  model.  First, 
for  sufficiently  low-energy  particles,  adding  more  points  becomes  impractical  because  efficiency 
gains  are  lost  and  the  cross  section  set  becomes  cumbersome.  Second,  even  a  large  number  of 
points  and  weights  does  not  help  to  resolve  the  angular  distribution.  Figure  5  presents  a 
comparison  of  angular  distributions  for  100  keV  electrons  in  gold  generated  using  a  32-point 
discrete  model  to  the  analog  model.  Discrete  artifacts  are  clearly  seen  in  the  figure.  For  this 
case,  adding  more  points  to  resolve  the  angular  distribution  will  result  in  a  MFP  that  is  on  the 
order  of  the  analog  MFP. 
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Figure  5.  Comparison  of  Angular  Distribution  for  1  MeV  Electrons  in  Gold. 

Regardless  of  incident  article  energy  and  atom  number,  the  GBFP  discrete  model  generates 
inaccurate  angular  distributions  when  using  lower  order  cross-sections  or  for  sufficiently  low- 
energy  particles  in  thin  slabs.  Therefore,  if  quantities  like  angular  distribution  or  energy 
spectrum  are  required,  a  hybrid  method  must  be  used  to  maintain  efficiency  while  mitigating  the 
discrete  artifacts.  The  discrete  artifacts  arise  because  particles  can  only  scatter  into  specific 
directions.  The  angular  distributions  tend  to  show  a  preference  for  these  particular  angles.  Use 
of  a  continuous  or  hybrid  model  will  resolve  the  discrete  artifacts.  Since  the  continuous  and 
hybrid  models  are  smooth  or  have  a  smooth  component,  the  discrete  artifacts  are  mitigated  [10]. 
This  is  the  motivation  for  using  the  hybrid  or  continuous  models.  However,  discrete  artifacts  are 
not  an  issue  for  calculation  of  integral  quantities  like  dose,  so  the  discrete  model  is  a  good 
candidate  for  dose  calculations. 
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The  table  below  summarizes  the  computer  time  required  to  complete  various  dose  calculations. 
As  indicated  in  fig.  1  and  in  table  3,  the  greatest  gains  are  achieved  for  higher  energy  source 
particles.  However,  the  predicted  speed-ups  from  fig.  1,  which  are  estimated  using  MFP,  are 
roughly  a  factor  of  10  higher  than  the  actual  speed-ups  based  on  computer  time. 


Table  3:  GBFP  Efficiency  Gains  Over  Analog  Monte  Carlo. 


Material 

Energy  (MeV) 

Slab  Width  (mm) 

Analog  Calculation 
Computer  Time 
(min) 

GBFP  Speed-up 
(Analog/GBFP) 

Gold 

1 

0.1 

147.1167 

19.1061 

Gold 

0.1 

0.1 

16.3000 

3.4077 

Aluminum 

1 

5 

591.0833 

51.3986 

Aluminum 

0.1 

0.5 

76.5500 

14.8161 

Au/Al 

0.25 

0.004/0.096 

54.1667 

9.3660 

It  is  important  to  note  that  each  dose  calculation  was  completed  using  only  two  discrete  angles 
and  one  discrete  energy  loss,  yet  there  is  still  good  agreement  between  the  analog  and  discrete 
solutions,  as  shown  in  Figures  6  through  13.  Further  accuracy  can  be  gained  by  inclusion  of 
more  points  and  weights.  Future  reports  will  include  details  on  the  impact  of  the  inclusion  of 
more  points  and  weights  on  the  accuracy  and  efficiency  of  dose  calculations. 

In  addition,  dose  profiles  were  calculated  for  a  target  with  a  high/low  Z  interface  to  demonstrate 
the  accuracy  of  the  method  at  interfaces,  shown  in  Figures  14  and  15.  The  GBFP  method  is  a 
single-event  Monte  Carlo  method  so  boundary  crossings  algorithms  are  not  necessary.  A  single¬ 
event  method  implies  that  it  is  a  true  Markov  process.  Therefore,  what  occurs  at  a  collision  site 
does  not  depend  on  how  the  particle  arrived  at  the  site.  Moreover,  a  Markov  process  can  be 
interrupted  at  any  point  and  restarted  without  introducing  error  or  correlation.  For  this  reason, 
boundary  crossings  require  stopping  the  particle  at  the  boundary  and  then  resampling  distance  to 
collision  in  the  new  material.  This  is  not  true  of  the  condensed  history  algorithm  because  this 
approach  is  not  a  Markov  process.  As  a  result,  condensed  history  algorithms  require  specific 
algorithms  for  handling  boundary  crossings  and  in  some  cases  condensed  history  codes  will 
switch  over  to  analog  Monte  Carlo.  Once  the  particle  is  far  enough  from  the  boundary, 
condensed  history  is  resumed.  None  of  these  additional  features  are  required  for  the  GBFP 
method  or  single-event  Monte  Carlo. 
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Relative  Difference  in  Dose  (1- Model  Analog) 


Analog  Discrete  2  angles,  1  Energy-Loss 


Figure  6.  Comparison  of  Dose  Profiles  of  100  keV  Electrons  in  Gold. 


Figure  7.  Relative  Difference  and  Uncertainty  in  Dose  Profile. 
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Relative  Difference  in  Dose  (1- Model  Analog) 


Depth  (cm) 

Figure  8.  Comparison  of  Dose  Profiles  of  100  keV  Electrons  in  Aluminum. 


Figure  9.  Relative  Difference  and  Uncertainty  in  Dose  Profile. 
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Relative  Difference  in  Dose  (1- Model  Analog) 


Analog  Discrete  2  angles,  1  Energy-Loss 


Depth  (cm) 


Figure  10.  Comparison  of  Dose  Profiles  of  1  MeV  Electrons  in  Gold. 


Figure  11.  Relative  Difference  and  Uncertainty  in  Dose  Profile. 
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Relative  Difference  in  Dose  (1- Model  Analog) 


Analog  Discrete  2  angles,  1  Energy-Loss 


Figure  12.  Comparison  of  Dose  Profiles  of  1  MeV  Electrons  in  Aluminum. 


Figure  13.  Relative  Difference  and  Uncertainty  in  Dose. 
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Analog  Discrete  2  angles,  1  Energy-Loss 


Figure  14.  Dose  Profiles  of  250  keV  Electrons  in  Slab  with  Gold-Aluminum  Interface. 


Depth  (cm) 


Figure  15.  Relative  Difference  and  Uncertainty  in  Dose  Profile. 
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5.  CONCLUSIONS 


Significant  progress  was  made  towards  implementation  of  the  GBFP  method  in  Geant4.  Both 
the  angular  distribution  results  and  the  dose  calculations  indicate  that  the  implementation  was 
done  correctly  and  that  Geant4  is  suitable  for  continuing  GBFP  research.  In  the  case  of  angular 
distributions,  it  was  demonstrated  that  the  approximate  solution  approaches  the  analog  solution 
by  addition  of  more  points  and  weights  when  approximating  the  cross-section.  This  behavior 
was  expected  and  indicates  that  the  algorithm  for  generating  the  cross-section  data  and  the 
implementation  of  the  model  in  Geant4  was  successful.  The  agreement  between  the  analog  dose 
calculations  and  the  discrete  dose  calculations  provide  further  evidence  a  successful 
implementation  of  the  discrete  model. 

Geant4  contains  the  necessary  infrastructure  for  modeling  a  wide  variety  of  problems.  The 
ability  to  model  sophisticated  geometries  will  enable  us  to  extend  this  work  to  multi-dimension 
problems  that  have  not  been  explored.  In  addition,  the  existing  multiple  scattering  and  ionization 
processes  and  models  will  be  used  to  compare  the  overall  performance  of  the  GBFP  method  to 
methods  within  Geant4  that  are  similar  to  condensed  history. 

Currently,  analog  and  discrete  GBFP  models  are  implemented  in  Geant4.  Implementation  of  the 
continuous  and  hybrid  models  will  begin  over  the  next  reporting  period.  In  addition,  secondary 
electron  production  is  not  accounted  for  under  the  current  GBFP  model.  Additional  work  must 
be  completed  to  introduce  secondary  production  and  correlation  of  energy-loss  and  angular 
deflection  for  inelastic  collisions  with  atomic  electrons.  Once  each  model  is  implemented  and 
tested  a  consistent  process  for  comparing  geant4  standard  multiple  scattering  and  ionization 
models  will  be  established.  In  addition,  simulation  of  sophisticated  geometries,  like  the  CEASE 
detector,  will  begin  once  the  models  are  implemented  and  tested. 
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